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Abstract 

We analyze non-integrable deformations of two-dimensional N = 1 supersymmetric 
quantum field theories with kink excitations. As example, we consider the multi- 
frequency Super Sine Gordon model. At weak coupling, this model is robust with 
respect to kink confinement phenomena, in contrast to the purely bosonic case. If 
we vary the coupling, the model presents a sequence of phase transitions, where 
pairs of kinks disappear from the spectrum. The phase transitions fall into two 
classes: the first presents the critical behaviors of the Tricritical Ising model, the 
second instead those of the gaussian model. In the first case, close to the criti- 
cal point, the model has metastable vacua, with a spontaneously supersymmetry 
breaking. When the life-time of the metastable vacua is sufficiently long, the role 
of goldstino is given by the massless Majorana fermion of the Ising model. On the 
contrary, supersymmetry remains exact in the phase transition of the second type. 



1 Introduction 



In two dimensions there are several bosonic field theories that have kinks as the basic 
excitations of their spectrum, the simplest example being the Ising model in its low 
temperature phase (see, for instance, pQ). Other examples are provided, for instance, 
by the minimal models Aim {m = 3,4, . . .) of Conformal Field Theory perturbed by the 
operator $13 [2]: these are off-critical theories with [m — 1) degenerate vacua, connected 
by massive kinks. 



Figure 1: Effective potential U{(p) of a quantum field theory with kink excitations. 

At the lagrangian level, theories with kink excitations are described by a scalar real field 
(p{x) with an action A = j (Px C and a Lagrangian density 
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^ = -f/M , (1.1) 
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where the potential U{lp) possesses a set of degenerate minima at vj^^"* {k = 1,2, ... ,n 



with ^^^li < V5^°'' < V^i+i- While these minima correspond to t 

I 0)fc of the associate quantum field theory, the kink stated | K^i (0) ) can be put in 
correspondence with the static solutions of the classical equation of motion 

dlcfix) = U'Mx)] , (1.2) 

with boundary conditions (p{—oo) = (pf^ and (p{+oo) = ^pf\ where I = kiil. The semi- 
classical quantization of these solutions has been discussed in [Sj HJ [5] . Conventionally 
I -^fc,A;+i(^) ) denotes the kink state between the pair of vacua | 0)^, | Q)k+i while the 
corresponding anti-kink is associated to | Kk+i,k{0) ). 



le different vacuum states 



^The quantity 6 in the symbol | Kki{9)) is the rapidity variable of the kinks, parameterising their 
relativistic dispersion relation E = M^-i cosh0, P = Mf^i sinh6', where Mf^i is their mass. 
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A popular example of such a lagrangian theory is the Sine-Gordon model, described 
by the potential 

U{ip) = —fi cos ay? , (1.3) 

that has an infinite number of degenerate vacua at ifn^ = 2TTn/a. The classical kink 
configurations that interpolate between two nearest ones are given by 



^±{^) = — arctan[exp(± a i//!^)] 



a 



(1.4) 

oo to X = +00, 



(+ refers to the kink whereas — to the antikink), so that, going from x = 
the field has a jump ±27r/a. 

The kinks obviously owe their existencies to the n-fold degenerate vacuum structure 
of the theory. However, such a degeneracy is often a fragile condition, which can be 
easily broken by inserting additional operators into the action. For instance, by adding 
an extra trigonometric interaction A cos unp (with u ^ a) to the Sine-Gordon model, the 
landscape of its potential changes from the situation (a) to the situation (b) of Figure 
2, even for small value of the coupling constant A. The unbalance of the vacua has a 
drastic consequence on the particle content of the theory: once the degeneracy of the 
original minima is lifted, some kink excitations (if not all) disappear from the spectrum 
of the perturbed theory, i.e. they get confined. As a consequence, the linear confinement 
potential between the kink and antikink causes the collapse of this pair into a string of 
bound states P [3 El [9l [10] (see Figure 3). 
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(a) (b) 

Figure 2: Potential of (a) Sine- Gordon model; (b) Perturbed Sine- Gordon model. 

Apart some special cases (some of them discussed in Section [22]), the above situation 
seems to be the general scenario for the bosonic theories with kink excitations. The situ- 
ation changes, however, if we add fermions. In fact, as shown below, adding fermions to 
a bosonic theory, while requiring a supersymmetry of the corresponding action, has the 
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effect to stabilize tlie vacuum states. In otfier words, supersymmetric tfieories are less 
sensitive to the confinement phenomena of the kinks, if put under a weak coupling defor- 
mation of the original action. At strong coupling, instead, they may present interesting 
phenomena of metastable vacuum states, which signal certain phase transitions. The two 
situations can be related to the spontaneous breaking of the supersymmetry that occurs 
at certain vacua. 

The paper is organised as follows. In the next section we briefly summarise the 
confinement phenomena of kinks in pure bosonic theories. In section 3 we introduce 
the supersymmetric theories and discuss their conformal limit. Section 4 deals with the 
Super Sine Gordon model, i.e. an integrable quantum field with kink excitations. In 
Section 5 we address the stability of the kinks under the supersymmetric deformation of 
the action that leads to the multi-frequency Sine Gordon model. In Section 6 we discuss 
the phenomena of meta-stable states that occur at finite value of the coupling constant 
of the perturbing operator. Our conclusions are presented in Section 7. 




(a) (b) 

Figure 3: (a) Kink- antikink state at distance x, where AE is the gap of the unbalanced 
vacua. The tendency of the system to shrink, in order to minimise the energy of this 
configuration, gives rise to the linear potential (b), with slope given by AE. The dashed 
lines are the bound states of the kink-antikink pair. 

2 Kink Confinement in Bosonic Theories 

Quantum field theories with kink excitations arc not necessarily integrable: (^^ theory 
in its Z2 broken phase, for instance, is not an integrable model although it has kink 
excitations. At the same time, the confinement of the kinks is also a general phenomenon, 
i.e. it occurs in deformation of both integrable and non-integrable theories. In order to 
have a certain analytic control of this phenomenon, it is however more convenient to 
discuss the confinement of the kink excitations that takes place in perturbed integrable 
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theories. The reason is that, for these theories, we can rely on a well-defined perturbation 
approach, the so called Form Factor Perturbation Theory (FFPT) 13, [H], that allows us to 
analitically follow the fate of the kink states: properly interpreted, the results obtained 
in the context of perturbed integrable theories turn out to be useful also to understand 
the confinement phenomena of a generic model. 

Based on what stated above, let us consider an integrable theory with purely bosonic 
degrees of freedom described by the action ^o- We assume that such a theory has n de- 
generate vacua I 0)fc, {k = 1, 2, . . . , n) and therefore its spectrum contains the topological 
excitations | Kk^k±i{d)) (the theory may also have bound states thereof). All kinks have 
the same mass M. 

Suppose now that we perturb the system using a field T(x): does the perturbed action 

^ = ^0 + A j (fx T{x) , (2.1) 

still support the topological kinks | Kk,k±i{6) ) as asymptotic states? As shown in [7l|8], 
to answer this question is sufficient to compute the correction 5Mki to the mass of the 
kinks due to the field $(a;). At the first order in A, this correction is given by 

5Ml ^ 2AFj(m) , (2.2) 

where 

F^i{9) ^ kmmKui{er)Kiu{e2)) (2.3) 

is the kink-antikink Form Factor of the operator T(a;), with 6 = 61—62. In integrable field 
theories (and this is the important technical point for which it is useful to consider these 
theories). Form Factors of a generic scalar operator 0{x) can be computed exactly. In 
fact, they satisfy manageable functional equations in virtue of the simple form assumed 
by the unitarity and crossing symmetry equations [121 fT3] . 

Consider the two-particle case: shortening the notation and denoting the kink state 
I -^A;,fc+i) by a and the anti-kink state | -ft'fc+i^fc) by a, for the Form Factor of the scalar 
operator O we have the following equations 

FgiO) = St{e)FP{-6) , (2.4) 

Fg{6 + 2z7r) = e-2-^°- FZ{-6) . (2.5) 

In the first of these equations the sum on h is on all kink states which start from the 
vacuum | 0)^, i.e. it includes the two states | Kfc.fcii)- It expresses the fact that in 
an integrable theory the two-particle threshold is the only unitarity branch point in the 
plane of the Mandelstam variable s = {pa + PaY = 4M^ cosh^(6'/2), the discontinuity 
across the cut being determined by the two-body scattering amplitude S^^^. 
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In the second equation the exphcit phase factor e '^™"io,a jg inserted to take into 



a 



When 



account a possible semi-locahty of the kink with respect to the operator 0{x) 
1o,a = 0, there is no crossing symmetric counterpart to the unitarity cut but when 
70,0 7^ 0, there is a non-locahty discontinuity in the s-plane with s = as branch point. 
However, in the rapidity parameterisation, there is no cut because the different Riemann 
sheets of the s-plane are mapped onto different sections of the 6'-plane; the branch point 
s = is mapped onto the points 6 = ±i7r which become, in this way, the locations of 
simple annihilation poles, whose residues are given by [I3] (see also [T6] ) 



-tResg=±,^Fg{e) = (l-e^2..7o,.) ^{o\0\0), . (2.6) 

Notice that these residues are expressed by the product of two terms: the first, determined 
by the semi-local index, while the second given by the vacuum expectation value of the 
field 0{x). 

With the above information, let's come back to the perturbed action (12.11) . It is easy 
to see that there is a drastic effect on the kink spectrum if (a) the perturbing operator 
T(x) is semi- local with respect to them and (b) if its expectation value of the vacuum 
I )a: is different from zero. 

If both conditions are fulfilled, the application of eq. (12. 2 p leads to an infinite correction 
to the kink masses. This divergence has to be interpreted, in turn, as the technical 
manifestation of the fact that the kink starting from the vacuum | )fc no longer survive 
as asymptotic particles of the perturbed theory. In order to clarify better this point, it 
may be useful to discuss a couple of example; 



□ry. 

I. 



2.1 Low-temperature Ising model in a magnetic field 

The simplest example of the above scenario is given by the Ising model. In its low- 
temperature phase it has two degenerate vacua, connected by a kink and an anti-kink 
excitation. Suppose now this system is coupled to an external magnetic field h, so that 
the perturbed action is given by {Ao is the integrable action of the low-temperature Ising) 

A = Ao + h j £xa{x) . (2.7) 

The magnetization operator a{x) is non-local with respect to the kinks, with non-local 
index given by 7^,0 = | [HI [IS] . Hence, the two-particle Form Factor of the operator 
cr(x) should have a pole at ^ = iivr. This is confirmed by the exact expression of this 

^Consistency of eq. (12. 5p requires 70^5 = — 7o,a- 
^It may be also useful to see [l4] . 
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matrix element, given (up to normalization) by [16j 

(0 I a(0) I a(5i)a(^2)) = tanh ^ ^ ^ . (2.8) 

Consequently, the kinks disappear from the spectrum of the perturbed theory (12.71) . 
producing though an infinite tower of bound state£| [SI [3, El E], [TD]. The confinement 
of the kink is evident from Figure 4: no matter how small the magnetic field h may 
be, the two original vacua are no longer degenerate in the perturbed theory. Hence, the 
perturbed theory (12. 7p does not have any longer kink excitations. 




2.2 Multi-frequency Sine-Gordon model 

As a second example, let's consider the perturbed Sine-Gordon model 

^ = ^(c^^v')^ + /i cosav^ + A cos(/5v9 + 5) . (2.9) 

This system has been analysed in detail in [S] , so here we briefly remind some of the main 
results. First of all, the quantum model does not depend separately by /i and A, it rather 
depends on the dimensionless variable 

In terms of ?7, we can identify two perturbative regimes of the model (12.91) where it is 
possible to apply the FFPT: the first regime is obtained in the limit ?7 — ^ 0, while the 

^ Not all of them are stable: the number of the stable ones changes by moving in the plane with axes 
(h,T) [nHH]: near the magnetic axis, the decay processes of the unstable particles can be computed by 
using the Form Factor Perturbation Theory [llj . 
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other is reached for — > cx). In the last case, one has simply to swap the role played by 
the two operators. 

Let us consider, for instance, the system in the vicinity of rj = 0. In the unperturbed 
situation {rj = 0), the field undergoes to the following discontinuity ip — > {p + 2TT/a across 
a soliton configuration. As a consequence, the exponential operator e*^'^*^^-' is semi-local 
with respect to the kink with semi-locality index 7^,5 = P/a. According to eg. (12.61) . 
the kink-antikink Form Factor of the perturbing operator T = cos{(3(f + 6) contains 
annihilation poles at ^ = iivr, with residues given by§| 

-iRese=±inFaa{0) = [cos6 -cos{6T'2Trf3/a)] (0|e^'^'^|0) . (2.10) 

Hence, for generic values of the frequency (3 (and the phase-shift 6) of the perturbing 
operator T(x), the above residues are different from zero. Consequently, the kink and the 
anti-kink of the original Sine-Gordon model become unstable excitations of the perturbed 
theory. This result that can be easily understood by noting that the perturbed lagrangian 
(12. 9p loses its original 27r/a-periodicity (see Figure 2). 

The possibility to control the changes of the spectrum in both perturbative limits 
{1] and 1] 00) allows one to deduce interesting information about its evolution in 
the intermediate, non-perturbative region. In certain cases, for instance, there may be 
a topological excitation in one limit which is no longer present in the other. When this 
happens, the very nature of topological excitations requires that a change in the vacuum 
structure of the theory takes place somewhere in the non-perturbative region, namely 
that a quantum phase transition occurs. Lines of phase transition are then expected 
to appear in the multi-frequency SG model for particular values of the parameters: as 
discussed in [8], this circumstance occurs when \6\ = n/n and P/a = 1/n. Once these 
values of the parameters are fixed, by varying t] a change of the vacuum state takes place 
at rjc'. at this critical value, the quantum field theory may be regarded as associated 
to a massless Renormalization Group flow that interpolates between the conformal field 
theory with central charge c = 1 (in the ultravioled regime) and the one associated to 
the Ising model, with c = 1/2 (in the infrared regime). The simplest way to see this, 
is to note that, at t]c, the effective theory at the vacuum nearby the origin is given by 
the massless 0^ associated to the class of universality of the Ising model: for rj slightly 
smaller than rjc, there are two degenerate ground states separated by a potential barrier, 
for 1] slightly greater than rjc there is only one vacuum with a parabolic shape, while at 
rj = rjc the quadratic term vanishes, giving rise to an effective 0^ behaviour with massless 
excitations (see Figure 5) 

^We choose as a vacuum state of the unperturbed theory the one at the origin, characterised by 
(0|(y5|0) ~ 0. With this choice the unperturbed theory is invariant under the reflection ip — > —ip and 
therefore it holds the equality (0|e*''^|0) = (0|e-^''^|0). 
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(a) 



(b) 



Figure 5: Potential of the multiple Sine-Gordon model with j3/a 
the vicinity of rjc-' (a) rj < rjc, (b) rj = rjc and (c) rj > rjc- 



(c) 

1/3 and S = 7t/3 in 



3 Bidimensional Supersymmetric Theories 

In this paper we will consider = 1 supersymmetric theories. They describe the dy- 
namics of a real superfield $(a;, 9) that consists of a scalar field (p{x), a real (Majorana) 
spinor ipa{x) and a real auxiliary field F{x), with the expansion 



^x,e) = ip{x) + eip{x) + leeF{x) 



(3.1) 



The space coordinates x^ — {x^,x^) and the the two real Grassmann coordinates 9^ — 
(^1,^2) describe the N — 1 superspace. We adopt the following representation for the 
7-matrices 



7° = 



The role of 75 is played by 



-i 

1 



75 = 7° 7^ 
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1 
-1 



i 

1 



(3.2) 



(3.3) 



The charge-conjugation is implemented by the matrix Ca/s — — (7°)a/3; so that for a 
charge-conjugate spinor tp^ = Ca/3'ipi3 — ip^, with ■0 = tp^^'^. Hence, in the above 
representation, the Majorana spinors and 9 are real. 
Under a translation in superspace 



x^^x'' + iCr9 , 9^^9^ + Ca 



(3.4) 



the superfield has a variation expressed by 



5^x,9) = LQa^x,9) , 



(3.5) 
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with Qa = d/dOa + i{l^d)ad^. The most general action invariant under the supersym- 
metric transformation (I3.5p is given by 



A = / d^xd^e 



where J d'^999 = 2, with the covariant derivative given by 



(3.6) 



^ i - (zd.Ye)^ . (3.7) 

ClUa 

Vr($) is the so-called superpotential, that we assume to be an analytic function of $. 
Integrating on the Grassmann variables, one arrives to the following expression of the 
action 



where W' {(f) = dW{ip)/d(p, etc. Finally, elimitating the auxiliary field F from its alge- 
braic equation of motion, i.e. substituting F —W'{(p) in the above expression, yields 
a lagrangian density for a supersymmetric theory given, in its most general form, by 

^ = 1 [(9,^)' + ^t'^M - llW'iv)? - lw'\^)^i; . (3.9) 
Associated to the transformation ( 13. 4p is the conserved supercurrent 

JH = (5.^)(7V^)a-^i^(7''^)a , (3.10) 

and the conserved supercharges 

Qa = j dx^Jl . (3.11) 
Let's also define the stress-energy tensor 

Tf"" = -i!-i^'^''iJ + ^^'^^''^--g^"'[{^o,^f -F"^] , (3.12) 
and the topological current 

= -e^-Fd,^ = e^^'d.Wiip) . (3.13) 
As shown by Witten and Olive [17J, the most general supersymmetry algebra is then 

{Qa, Q/S} = 2{^x)at3 + Z , (3.14) 

where = J dx^T^^ is the total energy and momentum, whereas Z is the topological 
charge 

2 = j dx'e = [W{v)tZ = Ui=+oo -W{v) Ui=_oo . (3.15) 



A convenient and explicit form way of the supersymmetry algebra (I3.14p is the following 



Ql = P+ , = P , Q+Q- + Q-Q+ = 2Z , (3.16) 

where we have used the notation (Qi, Q2) = {Q~, Q+) and P± = Po± Pi. 

In the above formulas, Z is different from zero only when admits several so- 

lutions of the equation W' {^) = 0, i.e. if the model supports kink excitations that 
interpolate between different vacua at = ±00. This is also evident from the expression 
of the topological current this is a purely derivative term which is different from zero 
only on field configurations that do not vanish at infinity. When it is evaluated on the 
classical expression of the kinks, Z is equal to their mass (see eq. fl4.32p below). Some 
examples of theories with non-zero Z will be considered in Section |H 

Eq. (13. 9 p gives the most general lagrangian of a supersymmetric theory. As evident 
from its expression, its interactions consists of two terms: the first one, a Yukawa coupling 
for the fermion {Vy = iy"(v^)^V'); the second one, a potential that involves only the 
bosonic degree of freedom (Vg = [ly '(</?)] ^). The last term is intrinsicaly positive and, as 
well-known [18], it can be used to study the spontaneously breaking of supersymmetry: if 
[Vr' ((/?)] 2 has zeros, they are the true vacua of the theory and supersymmetry is unbroken; 
viceversa, if [W'{ip)]^has local minima that are not zeros of this function, supersymmetry 
is spontaneously broken at these minima. They are called meta-stable vacua if exists a 
true vacuum somewhere in the landscape of (</)). The possibility of constructing reliable 
theories with a broken supersymmetry at these meta-stable vacua is a topic currently 
under an active investigation, after the seminal work [19]. We will further comment on 
this topic in Section [61 



3.1 Superconformal minimal models and deformations 

As for pure bosonic theories, supersymmetric models may present a critical behavior. In 
the following, our attention will focus on the so-called = 1 discrete series SMm [20] . 
identified by combining the conformal invariance of the critical point together with the 
supersymmetry and the unitarity of the model^. The central charge of these conformal 
theories assumes the discrete values 

c=|- , m = 3,4,5,... (3.17) 

2 mini + 2) 



^In this section we will only remind the properties of these models that are useful in the sequel of 
the paper. For a detailed discussion of the super-conformal models, in addition to [20], see for instance 
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The fields Or^s;f,s, with l<r, r<m — 1 and 1 < s, s < m + 1, have anomalous dimensions 
1] = Ar^s + Af^s, where 

\(m + 2)r — ms]'^ — A 1 r , -.^ d , 

Ars = A.m-rm+2^s = ^ + — 1 - -1 • 3.18 

' + 8m(m + 2) 32 ^ ^ ^ J ^ ^ 

When r — s is an even number, the field Or,s;f,s is a primary superfield that belongs 
to the Neveu-Schwartz representation of the super-algebra. When r — s is instead an 
odd number, the corresponding field Or,s;f,s is an ordinary conformal primary field that 
creates a representation of the Ramond sector of the super-algebra. 

As shown by Zamolodchikov [23], the N = 1 superconformal model SAim (with the 
operator content fixed by the diagonal partition function) can be identified with the 
Landau-Ginzbrug lagrangian of a real scalar superfield $ with the superpotential given, 
in this case, by W{^) = ^ 

The first superconformal model, with m = 3 and c = ^, corresponds to the class of 
universality of the Tricritical Ising Model (TIM) |20]. According to the above identifica- 
tion, this model is associated to the superpotential IV($) = | after eliminating the 
auxiliary field F, it gives rise to the interaction terms 

V = ]-^^ + ^^^ . (3.19) 

The next model, with m = 4 and c = 1, corresponds to the class of universality of 
the critical Ashkin- Teller model (i.e. a particular gaussian model). Its superpotential 
Vr(0) = 1$^ produces the interaction terms 

^ = ^^' + ^^^^' ■ (3.20) 

An important difference exists between the superconformal minimal models SM-m when 
m is an even or an odd number. To see this, it is useful to computed the Witten index 
Tr(— 1)-^ of the superconformal models (/ is the fermion number, which assumes the 
value / = for bosonic states and / = 1 for fermionic states). As usual, this can be 
computed by initially defining the theory on a cylinder [26]: the Hamiltonian on the 
cylinder is given hj H = = Lo ~ c/24, where Q is the holomorphic compoment of the 
supercharge, c/24 is the Casimir energy on the cylinder and Lq = ^^dzzT{z), with 
T{z) the holomorphic component of the stress-energy tensor. Since for any conformal 
state I a) with A > c/24 there is the companion state Q | a) of opposite fermionic parity, 
their contributions cancel each other in Tr(— 1)-^ and therefore only the ground states 
with A = c/24 (which are not necessarily paired) enter the final expression of the Witten 
index. For the minimal models, by using the conformal dimensions given in eq. fl3.18p . it 
is easy to check that there is a non-zero Witten index only for m even: the Ramond field 
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of these models has indeed conformal dimension A = c/24. Therefore the lowest 
superconformal minimal model with a non-zero Witten index is the one with m = 4, i.e. 
the gaussian model, whereas the Tricritical Ising Model has a zero Witten index. 

It is well known that the Witten index is related to the supersymmetry breaking 
problem: if it is nonzero, the supersymmetry cannot be broken and if the supersymmetry 
is broken it must be zero. This is explicitly confirmed by the examples studied in the 
literature. In the case of the Tricritical Ising Model, perturbing the critical action by 
the super-operator O13, i.e. inserting the extra term A $ in the supersymmetrix action, 
produces the bosonic interaction term (see Figure 6) 

Vb = li^' + X? . (3.21) 

If A < 0, the ground state energy is zero and supersymmetry is unbroken: both scalar and 
fermion fields are massive and the theory has kink excitations [27|. Viceversa, if A > 
the ground state energy is non-zero and supersymmetry is spontaneously broken: the 
scalar field if becomes massive whereas the fermion stays massless and plays the role of 
goldstino. In this case, the Landau- Ginzbrug theory describes the massless flow from the 
Tricritical Ising model to the critical Ising model [26]. This flow gives rise, in particular, 
to an integrable field theory with the exact ^-matrix determined in [28]. Notice that, at 
the critical point, Vb has a behaviour ~ y?'^, that obviously does not correspond to 
the class of universality of the Ising model: in fact, one has to take into account that 
the superconformal model SMls has also additional interaction relative to its fermionic 
degree of freedom, i.e. the nteractions of the model are not only expressed by Vb- 




(a) (b) (c) 

Figure 6: Shape of the bosonic interaction term of the Tricritical ising Model by varying 
\: (a) X < (susy exact) (b) X = (superconformal point) and (c) X > (broken susy) . 

In the gaussian model, the role of the super-operator (9i 3 is played by the term 
Perturbing the conformal action by means of this operator, one has (see Figure 7) 

Vb = ^^'((^' + A)2 . (3.22) 
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By varying A, the landscape of the bosonic potential changes as in Figure 7: in this case, 
no matter what the sign of A may be, the potential has always at least one zero, i.e. 
supersjTiimetry is always unbroken. Notice that at the critical point, Vb has a behaviour 
at the origin Vb ~ ^p^- For A < 0, the excitations of the model are massive kinks while, 
for A > 0, they consists of massive scalar and fermion particles localised around the exact 
vacuum of the theory. 



\r\/\l 



(a) 



(b) 



(c) 



Figure 7: Shape of the bosonic interaction term of the gaussian model by varying X: (a) 
A < (three degenerate vacua of zero energy: susy exact) (b) X = (superconformal 
point) and (c) X> Q (unique vacuum of zero energy: susy exact) . 



4 Kink in supersymmetric theories 

In the following, as an example of a supersymmetric integrable theory with kink excita- 
tions, we will consider the super Sine-Gordon mode (SSG). This model has been studied 
by several authors and a great deal is known about it: its quantum integrability, for 
instance, has been checked in [29], its exact S-matrix has been proposed in [30j and it 
follows the general classification given in [31], whereas its analysis in the framework of 
perturbed conformal theories has been carried out in [33] . Its thermodynamics and finite 
size spectrum has been studied in [32]. The long-standing and subtle issue of quantum 
correction to the mass of the solitons has been solved in the series of papers [HH [35] , 
where it was shown that the soliton states satisfy a BPS property (i.e. their mass being 
equal to their topological charge) even beyond the tree level. 
As action of the SSG we will take 



A = / d^xd^e 



i(D,<l>)D„$-4cos(«*) 
4 



(4.23) 



Eliminating the auxiliary field F, it gives rise to the Lagrangian 

^ = 1 [id.^y + ^^rd^^] - ^ sin2(a cp) - ^ cos(a if)^lJ^P . (4.24) 
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In the standard perturbative approach of this model, the parameter fi (that we assume 
to be a positive quantity) plays the role of the mass of the scalar excitations nearby each 
minima of Vb = ^sin^(av9), that are located at ipf^ = kn/a. All these minima have 
zero energy, i.e. they are all possible supersymmetric vacua of the quantum theory. A 
more detailed discussion on the nature of these vacuum states will be done at the end of 
this section. Notice that the above theory may be also regarded as a massive deformation 
of the free conformal action with central charge c = 3/2, given by a massless bosonic field 
and a massless fermionic fielcjl|. 

The theory ( 14. 24^ is invariant under the Z2 parity ip —ip and under the shift 

ip ^ ip + n — . (4.25) 

a 

It is also invariant under a half-period shift of the bosonic field 

^-^^ + n- , (4.26) 
a 

a pact to change the relative sign of the fermion by 75, ip —>■ 'j^ip, i.e. 

^i^^i , ^2^ -^2 . (4.27) 

In view of the periodicity fl4.25p . we can always focus our attention to the interval 
(0,27r/a) of the field <p. 

The minima of Vb are connected by (supersymmetric) kinks | S), which can be 
identified by the conditions [IT] 

(Q+±Q_) I 5) = , (4.28) 

where the + refers to the kink whereas the — to the anti-kink. For the bosonic component 
of these kink states, the above conditions end up in solving the first order differential 
equations^ 

^ = ±F(^'') = ±W'(ip'') . (4.29) 
dx 

The kink and the anti-kink solutions (corresponding to + and — respectively) of the SSG 
model are explicitly given by 

(^5|(a;) = - arctan(e±''^) . (4.30) 
a 

They connect the adjacent vacua (p^^^ = and (p>^^^ = n/a, all other kinks or antikinks of 
the model being equivalent to the above. Their classical mass M is expressed in terms 
of their topological charge, thanks to the identity 

p+ + p_ = (g+ ± g_)2 T 2z (4.31) 



''The discussion of this point of view can be found in [33j . 

*See for an elegant discussion of the classical kink solutions in terms of the superspace formalism. 
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that follows from the supersymmetry algebra (I3.16p : for the kink (anitkink) at rest, the 
left hand side is given by P+ + P_ = 2M, whereas the right hand side, using eq. (14.281) . 
is equal to T'^Z, therefore [IT] 

M = \Z\ . (4.32) 

As a matter of fact, the kink and the anti-kink are double degenerate states [37]. One 
way to see this is to observe that, in the presence of kink solutions, the Dirac equation 
for the fermion field 

il^d^i)- II cos{aifif) ip = (4.33) 

possesses static normalizable solutions \ localised at the position of the kinks [36l 1371 
[38] . Namely, they satisfy the static differential equation 

i7^5iV±^ -/xcos(a(/?^)^i°^ = . (4.34) 

For the kink backgrounds y?^ of the SSG model, their explicit expression is given by 

(A^ is a finite normalization constant). These are zero-energy eigenf unctions of the Dirac 
equation, as it can be seen by looking for solutions of eq. fl4.33p in terms of eigenfunc- 
tions, ip{x,t) = e~^'^''* ipplx) . Let's focus the attention, for simplicity, only on the kink 
configuration, since for the antikink one has an analogous situation. We can expand the 
fermionic field in terms of its eigenfunctions: taking into account that ip is a Majorana 
field, one has 

^ = ao^f + J2 (^P V + i^r,-) ■ (4.36) 
p 

The operators oq and bp satisfy anticommutation relations 

{ao,aJ} = 1 , {bp,bl} = 6p,g , (4.37) 

and other anticommutation relation vanish. The operators bj^ and bp are standard cre- 
ation and annihilation for the non-zero energy fermions in the presence of the kink. The 
operator oq, associated to the fermion zero mode, has instead a bidimensional represen- 
tation: operating on the kink state, it produce another kink state with the same energy 
but different fermion number. Hence, in supersymmetric theories, the kink is a doublet. 
Denoting these two states by | A' ; ± ) , one has 

ao\K;-) = , \ K ; +) = al \ K ; -) . (4.38) 

The same considerations also apply to the antikink state. 
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In view of the above discussion, the nature of the vacuum states of the SSG model is 
as follows: every vacuum y^g^^ = is non-degenerate, while those located at y^fj^i = 
(2n+i) TT double degenerate. As discussed in |33j on the basis of an argument by 
Zamolodchikov [27j, this can be also argued by the sign of the effective mass term of the 
fermion field 

W = I cos{aip)ijij (4.39) 

at these vacua: when the effective mass is positive, the vacuum is not degenerate, while 
is double degenerate if the effective mass is negative. Hence, altogether we have the 
situation shown in Figure 8. 




5 Mult i- frequency Super Sine-Gordon Model 



Let's now deform the supersymmetric action (14.231) by inserting another trigonometric 
term, i.e let's consider the action 



A 



(fx (fe 



-iD^^) D„$ - ^ cos(a<l>) - 4 cos(/3$) 
4 a p 



(5.1) 



where we have opportunately rescaled the couplings. By using the super-Coulomb gas 
formalism [31] and the operator product expansion of the (super) exponential operators 
present in (15.11) . one can establish the renormalizability (i.e. the stability) of this action 
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in the range of the parameters 



< Stt , /3^<87r , a (3 < Aw . (5.2) 

In the following we assume the validity of these inequalities. Moreover, we restrict the 
attention to the case in which the ratio uj of the two frequencies is a rational number 

uj = — = - , (5.3) 
a q 

with p and q two co-prime natural numbers, p < q. Since any irrational number can be 
approximated with arbitrary precision by the sequence of rational approximants given 
by its continued fraction expressions, the case when u is irrational may be regarded as a 
particular limit of the rational situation. 

Since the two trigonometric interactions enter the action (15.11) in a symmetric way, 
there are two natural perturbative regimes of the theorjj^: A — > (with /i fixed) or /i ^ 
(with A fixed). Below we will deal with the first regime, since the other can be simply 
recovered by swapping the role of the two operators. 

The question we would like to address concerns the fate of the kinks of the original 
SSG: are they going to be confined once the new interactions is switched on? Or, on the 
contrary, are they going to remain stable excitations of the perturbed action (15. ip ? To 
answer these questions, let us first of all eliminate the auxiliary field F and then focus 
the attention on the resulting bosonic interaction Vb associated to the action (15. ip 

Vb{^, ^) = sm{a^) + A sm{P^)f . (5.4) 

To simplify the following formulas, let us rescale the field (f ^ ip/a and the coupling 
constant A ^ A//i in such a way that, up to an overall constant, Vb assumes the form 



VB^f,^) = [sin(/) + A sin (c<j </))]^ 



sin 09 + A sin - 09 
q 



1 2 

p ■ 



(5.5) 



The first observation is that, no matter what are the values assumed by A, the origin 
(^(0) = is always a zero of this expression, i.e. this model always possesses a true 
supersymmetric vacuum. We will comment later on the possibility of having meta-stable 
vacua, i.e. vacua where the supersymmetry is spontaneously broken, by varying A. 

The second observation is that, at lowest order in A, all the vacua of the original SSG 
potential continue to remain degenerate, i.e. the original kinks are always stable at weak 
coupling! The simplest way to show this, is to follow the evolution of the zeros of (15. 5p 



^As in the bosonic case [5], the quantum theory actually depends on the appropriate dimensionless 
combination of the two couplings. 
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when we switch on A. At A = 0, they are placed at (fn'^ = nrc. Switching on A, we can 
look for the new location of the zeros in power-series in A, i.e. 

^ nvr + Aei^) + . . . (5.6) 

Substituting this expression into the equation Vb(v^, A) = 0, it is easy to prove that, 
at the first order in A, Vb(v^, A) has the same zeros of Vb(v^,0). This can be seen as 
a simple consequence of the fact that Vb is given, in supersymmetric theory, by the 
square of a function: namely, if Vb(x) = [/(s)]^ and f{x) has zeros at x = Xi,X2, ... in 
the interval /, perturbing this function by Xg{x) one has Vb{x,X) = [f{x) + Xglx)]"^ ~ 
[/(x)]^ + 2Xf{x)g{x). Therefore, if g{x) does not have zeros in the interval /, at the 
lowest order in A, Vb{x, A) has the same zeros of Vb{x) in /. The actual shift of the zeros 
can be computed by imposing that VB^f, A) has zeros also at the second order in A: the 

consistency of this request comes from the result there is always a solution in terms of 

(1) • u 
eil , given by 

ei^) = (-irsin!^ . (5.7) 

Hence, at weak coupling, the zeros (fn^ shift their position by en but they do not disap- 
pear. 

The stability of the kinks at the lowest order in A and the role played by supersymme- 
try becomes even more evident if we apply the considerations of Section [2l For doing so, 
the first thing is to identify the operator that, at lowest order in A, spoils the integrability 
of the original model. It is easy to see that the role of T is played in this case by the 
operator 

6 = sinip sin^ujip) + ujipilj cos{unp) . (5.8) 

To compute its semi-local index 7e „ on the kink state, we can use the fact that the 
fermionic and the bosonic field are components of the same superfield and that the 
exponential operators e*'''^ has a semi-local index 7 = ^ with respect to the kink of this 
theoryF^I. If we quantize the theory by choosing as vacuum the zero at the origin, for the 
residue of the Form Factor of 6 that controls the confinement of the kinks we have 

Rese=±in Fg{e) = [1 + cos(7r^)] o(0|e(0)|0)o . (5.9) 

Notice that the above quantity consistently vanishes if u; = 1, i.e. if we had decided to 
perturb the action by the same trigonometric term that enters its original definition: in 
this case the deformation just redefines the original coupling constant fi and this cannot 
have consequences on the confinement of the kink. For all other values of u, the term 



^"Notice that this is half the semi-local index of the bosonic Sine-Gordon case, since the supersymmetric 
bosonic potential sin^ ip has a frequency that is twice the frequency of the usual Sine-Gordon model. 
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[1 + cos{ttuj)] is different from zero but the residue is nevertheless zero, because it is the 
vacuum expectation value of G that vanishes in this case! This is a consequence of the 
unbroken supersymmetry of the theory: in fact, G is nothing else that an additional term 
in the trace of the stress-energy tensor of theory and, since the supersymmetry is exact, 
from Q-t I 0)0 = simply follows o(0|G(0)|0)o = [ID]. It is interesting to notice 
that the absence of the confinement of the kinks is also independent from the rational or 
irrational nature of u. 

From a naive point of view, the vanishing of the vacuum expectation value of the 
perturbing operator G is due to the simultaneous presence of a bosonic and a fermion term 
in its expression: in this case, the vacuum expectation value of the fermionic term tends 
to cancel the bosonic one. The exact cancellation obviously happens for an opportune 
fine-tuning of the relative couplings of these terms - a fact that simply expresses the 
supersymmetry invariance of the theory. In the usual multi-frequency Sine-Gordon, this 
cancellation is absent for the only presence of the bosonic term: the result is then, in 
general, the confinement of its kinks. 

The vanishing of the residue (15.91) implies the absence of confinement of the kinks at 
weak coupling but, obviously, this does not imply that their mass remain untouched. As 
a matter of fact, the masses change and, in particular, the kinks are no longer degenerate 
when A is switched on. Namely, the original kink excitations split into families of long- 
kinks and short kinks: the former overpass higher barriers, the latter lower barriers (see, 
for instance. Figure 9). The actual computation of the new masses of the kinks can be 
done by employing their topological charges, i.e. by using eqs. (13.151) and (14.321) . together 
with the shift of the zeros given by (15. 7p . 

6 Phase transition and meta-stable states 

The analysis done in the previous section applies in the weak coupling regime A ^ 0. One 
may wonder, however, what is the evolution of the spectrum!"] of the theory by varying 
A to finite values or even pushing it to the strong coupling limit A — 00. The result , as 
we will see, is quite interesting. 

Consider the bosonic and the fermionic potential terms of the theory, in their rescaled 
form 



^^The SSG model has, in addition to the kinks, also bound states thereof. In this paper, for the sake 
of simplicity, we focus our attention only on the kink excitations of the theory. 




(6.1) 
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Vy 



cos (p + X - cos 

p 



p 



In the presence of A, the periodicity of the theory is no longer equal to 27r but has become 
27rgf. Therefore it is convenient to plot Vb on this extended interval, as in Figure 9. 




(a) (b) 

Figure 9: Plot of Vb with A = 0.1 in the interval (— 27rg, 27rg) for q — 3: (a) p = 1 and 
(h) p = 2. 

The evolution of the spectrum by varying A is constrained by a simple argument. At 
A = 0, one has the sequence of the zeros of the first trigonometric term, which become 
slightly displaced from their original position when A is small, without changing though 
their number. Denoting by Nz{X) the number of zero at a given A in the enlarged interval 
(0, 27rgf), we have 

A^,(0) = 2g + 1 (6.2) 

(in the above number we have also included the zero at the origin). On the other hand, 
when A — > oo, the number of zeros becomes 

lim N,{X) = 2p+ 1 , (6.3) 

A— >oo 

i.e. the number of zeros in the interval (0, 27rg) of the second trigonometric term in Vb- 
Hence, by varying A, there should be a variation of the number of zeros equal to 

AN, = 2{q-p) . (6.4) 

Since the kinks own their existence to the zeros, a variation of their number implies that 
certain kinks will disappear by moving A from A = to A = oo. Since they are topological 
excitations, their disappearance signals that certain phase transitions will take place in 
the model at some critical values of the coupling, A = X'^'' . At these critical values, the 
system has massless excitations that will rule its long-distance behaviour. As discussed 
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.(1) 



.(2) 



N 



.(3) 



Figure 10: Nz{X) vs X, with the sequence of the critical values Ac 

below, A^"'' are identified as the values where A^^(A) jumps by a step of 2 units (Figure 
10). 

The way in which the number of zeros changes, by varying with continuity A, is 
that pairs of zeros collide and then move on immaginary values. When this happens, 
the barrier between the two colliding vacua vanishes and, correspondingly, the kinks that 
connect them become massless. Notice that, in the evolution of the zeros, the ones placed 
at (f^^^ = 0, yj^^-* = Tcq and (^^'^^ = 2Txq will never disapper. 

There is, however, a different way in which the disappearance of the zeros is imple- 
mented, according whether (p — q) is an even or an odd number. In fact, when (g — p) is 
an odd number, among the pairs of colliding zeros there will always be a couple of them 
placed just on the right and on the left position of (^9*^°^ = ixq: at some critical value Ac, 
these zeros will collide but strangling, in the meantime, the zero at vrg that is in between. 
When [q — p) is an even number, this will not happen. Hence, the conclusion is the 
foUowing^i: 

• When (g — p) is an even number, there will be a sequence of phase transitions that 
will recall the phase transition of the Tricritical Ising Model, i.e. a phase transition 
with a local spontaneously supersymmetry breaking (see eq. (I3.2ip and the relative 
discussion). At the critical points, the local form of the potential Vb at the position 
(p of the colliding zeros is given by ~ (<^ — (pY. 

• When (g — p) is an odd number, in addition to a sequence of TIM-like phase tran- 
sitions, at a certain value of the coupling, A = Ac, there is also a phase transition 
localized at the vacuum (f = irq. This phase transition will recall the phase tran- 
sition of the gaussian model, i.e. a phase transition without a local spontaneously 

^^The analysis refers to A > 0. If A is negative there is a swap of the two cases. 
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supersymmetry breaking (see eq. fl3.22l) and the relative discussion). At A = Ac the 
local form of the potential Vb at ip = irq is given by Vb ^ — vrg)^. 

In order to make more transparent the above analysis, it is useful to consider explicitly 
two examples: the first refers to = |, i.e. {q — p) = 2, the second instead to = |, 
with {q — p) = 1. 




(d) (e) (f) 

Figure 11: Plot o/Vb in the interval (0, Gtt) foru = 1/3 by varying X: (a) small value of 
A; (b) A < Ac; (c) A = Ac/ (d) A > (e) A > A^; A ^ oo . 

Let us consider first the case u = ^. By varying A, the landscape of Vb changes as 
shown in the sequence of the pictures of Figure 11. Observe that, just before reaching 
Ac, there are two zeros in the interval (0, Svr) that are going to merge together (obviously 
there are similar zeros in all other intervals obtained by periodicity). Let's call these 
vacua (p^^^ and (p^2^. At A = Ac the barrier between these two zeros vanishes and, soon 
after Ac, a meta-stable vacuum is created at their middle position ip = ipmv For A close 
to Ac, the shape of Vb at ipmv can be parameterized similarly to the bosonic potential of 
the Tricritical Ising Model 

Vb ~ [{if - Vn^vf + {\ - \c)Y . (6.5) 

Since at A = Ac some topological excitations disappear, a phase transition takes place in 
the model: if we had originally decided to quantize the theory either at (p^^^ or at (p'^\ the 
effective dynamics seen by any of these two vacua appears as a cascade of massless flows: 
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the first, at A = from c = 3/2 to c = 7/10 (tlie central cliarge of TIM) and soon after 
Ac, from c = 7/10 to c = 1/2, the last value being the central charge of the Ising model. 
This jumps in the central charge are consistent with the c-theorem [Hj and with the 
unitarity of the theory. Although the zero at the origin always leaves the supersymmetry 
of the model exact, the quantization around the meta-stable vacuum state at (pmv realises 
however a spontaneously supersymmetric breaking situation, with a spectrum given by 
a massive scalar and a massless fermion. This can be explicitly checked by studying Vy: 
for A > Ac (and for all the range of values of A that permits to consider the lifetime of 
this meta-stable vacuum sufficiently long), the term that multiplies ipip is indeed zero at 
= '^vm- The massless fermion is nothing else than the usual Majorana fermion of the 
critical Ising Model. Increasing further A and reaching values much larger than Ac, the 
lifetime of the metastable vacuum shortens and the effective theory based on this pseudo 
vacuum loses, at certain point, its validity. As a matter of fact, at strong coupling this 
pseudo-vacuum is going to be finally absorbed into the maximum of the second term of 
Vb, i.e. sin (Iv^). It is easy to check that the situation just discussed for the case uj = ^ 
also occurs in all other cases when (q—p) is an even number, possibly with several critical 
values of A where phase transitions of the type of the Tricritical Ising Model occur. 

Let us discuss now the second case, = §• The evolutation of the landscape of Vb by 
varying A is shown in the sequence of pictures of Figure 12. When A is just switched on, 
the number of zeros does not change but the barriers that separate them start to lower. 
If one concentrates the attention on the barriers placed in the middle of the plots, one 
notices that their evolution follows the patter presented by the gaussian model. Denoting 
by Ac the critical value of the coupling, in the vicinity of this value the shape of the Vb 
can be parameterized similarly to the bosonic potential of the gaussian model 

Vb ~ - 3nY [(^ - 3nf + (A - Ac) ] ' . (6.6) 

In this case, overpassing the critical value Ac, supersymmetry still remains exact at the 
vacuum (^^^^ = 3tt. It is easy to see that a similar situation occurs in all other case when 
[q — p) is an odd number, the only difference possibly being the existence of a sequence 
of other critical values Ac"'' (in addition to Ac) where TIM-like phase transitions occur at 
other vacua. This happens if {q — p) > 1. If we had decide to quantize the theory by 
choosing as vacuum the one placed at (p^^^ = vrg, the dynamics seen at this vacuum state 
consists, at A = Ac, in a massless fiow from c = 3/2toc=l, and then from c = 1 to 
c = 0, without a breaking of the supersymmetry. 
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(d) (e) (f) 

Figure 12: Vb foru = 2/3 by varying \: (a) small value of \ (b) A < Xd (c) A = Xd (d) 
X > Ac/ (e) X > Ac/ (f) X^oo . 

7 Conclusions 

In this paper we have analysed, through the study of a particular model, the breaking of 
integrability of a supersymmetric theory. On the contrary to purely bosonic theories, this 
does not lead to a confinement of the kink exitations. The absence of this phenomenon 
is due precisely to the unbroken supersymmetry invariance of the theory. However, by 
varying the coupling constant, the spectrum changes. 

At weak coupling, there is no longer the degeneracy of the original kinks: the kinks 
split into a sequence of long and short kinks. Since they are still BPS states, their mass 
can be computed in terms of their topological charge Z. An interesting open problem 
consists of comparing the values obtained by this approach with the masses computed in 
terms of the Form Factor Perturbation Theory. In order to do so, one needs to generalize 
the supersymmetric Form Factors analysed in P2] to theories with kinks. 

At finite values of the coupling, there exist critical points where phase transitions 
take place in the system: at these critical values, pairs of kinks become massless and 
then disappear from the set of asymptotic states. 

Among the phase transitions that take place in the theory, some of them realise a 
spontaneously supersymmetry breaking around certain meta-stable vacua. This happens 
if (q—p) > 1, where u = p/q is the frequency of the perturbing superfield. A quantization 
of the theory around the metastable vacua leads, nearby the critical values of the coupling, 
to the effective dynamics that describes the spontaneously supersymmetry breaking of 
the Tricritical Ising model. This consists of a massless flow from this model to the critical 
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Ising model, going from the short to the large distance scales. In this flow, the scalar 
particle is massive while the fermion is the massless Mayorana particle of the Ising model, 
i.e. the goldstino of the supersymmetry breaking. 

If (g — p) is an odd number, the system also present a phase transition at Ac which 
does not break supersymmetry at the vacuum in the middle of the interval (0, 27rq) 
where it occurs. In this case, the effective dynamics seen at this vacuum consists of the 
class of universality of the gaussian model, which presents massive and supersymmetric 
excitations either before or after the critical point, becoming massless just at A = Ac: for 
A < Ac they appear as kinks, for A > Ac, they are instead ordinary particles. 
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